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The second law of thermodynamics imphes that no macroscopic system may oscillate indefinitely 
without consuming energy. The question of the number of possible oscillations and the coherent 
quality of these oscillations remain unanswered. This paper proves the upper-bounds on the number 
and quality of such oscillations when the system in question is homogeneously driven and has a 
discrete network of states. In a closed system, the maximum number of oscillations is bounded by 
the number of states in the network. In open systems, the size of the network bounds the quality 
factor of oscillation. This work also explores how the quality factor of macrostate oscillations, 
such as would be observed in chemical reactions, are bounded by the smallest equivalent loop of 
the network, not the size of the entire system. The consequences of this limit are explored in the 
context of chemical clocks and limit cycles. 
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I. INTRODUCTION 

Oscillations are ubiquitous. We celebrate them and at- 
tempt to harness them. Naturally, this interest drives us 
to study them. There has been no shortage of analysis 
of the simple harmonic oscillator in all of its variations, 
but much of the periodicity around us is not equivalent 
to a mass on a spring. For example, the beating heart 
is driven by molecular motors which exist in the low- 
Reynolds number regime, where viscous damping over- 
whelms inertia. For these molecular constituents, buf- 
feting by solvent molecules prevents coherent oscillations 
from persisting on a timescale longer than the mean time 
between collisions, which is on order picoseconds [1]. De- 
spite this, we observe the coordination of overdamped 
components to produce periodic behavior [2 . Studying 
this coordination on a problem-by-problem basis has un- 
covered some conceptual principles for creating oscilla- 
tory behavior in the overdamped regime, but few truly 
fundamental laws exist ^-^7| . This work bounds the oscil- 
latory performance of all discrete-state overdamped sys- 
tems, providing a new look at the necessary conditions 
for creating coherent oscillations in dynamic networks [8 . 

When the energy landscape of an overdamped system 
can be divided into distinct basins of attraction with bar- 
riers higher than /c^T, the system will tend to reach 
a local equilibrium within a basin before fluctuations 
stochastically drive it over a barrier into a neighboring 
basin. When this condition holds, it is common and ap- 
propriate to model each basin as a distinct state, with a 
static rate of transitioning from one state to another [9 - 
[11], meaning the entire system may be represented as 
network with nodes representing each state and a edge 
connecting nodes with a non-zero transition rate. These 
systems are finite state first-order Markov processes [12] 
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and can be modeled by the finite-state master equation: 

^ = ^n^p^it)-Y^nMt), (1) 

where Tij is the transition rate from state j to state i. 
For introductions to the master equation and its numer- 
ous applications, see [TSHIS]. We will assume that all 
rates are time-independent, meaning no external factors 
change the rates (but this does not necessarily mean that 
the system is closed). We also make the assumption that 
T is an irreducible matrix, enforcing the trivial condition 
that we are not modeling multiple mutually isolated sys- 
tems. Finally, we assume that the systems conserve prob- 
ability, which can always be enforced by adding states to 
the system to represent sinks or sources. The solution to 
Eq. [l]is p{t) = exp(Tt)p(0), where T is matrix notation 
for Tij, Tii = — ^j-^^Tji and p{t) is vector notation for 
Pji'^) Hz]- Systems represented by the master equation 
are completely described by the transition rate matrix, 
T, and the initial conditions p{0). The complete solution 
is 

m = J2 Vje'^' • P(0) + a, (t)) , (2) 

j 

where vj is the j^^ eigenvector, is the inverse of 

the matrix of eigenvectors, and aj{t) is a polynomial due 
any eigenvalue degeneracy [18]. Hence, characterizing the 
properties of T also characterizes the dynamics of the sys- 
tem Q3l[T3|20|. Because the time dynamics of individual 
modes are ultimately determined by the eigenvalues of T 
(see Eq. [2|, we will be concerned with these eigenval- 
ues and how they relate to oscillations. We first explore 
these eigenvalues and prove how they constrain the pos- 
sible oscillations to be fewer than the number of states in 
the system. Second, the present work provides examples 
of these limits by exploring the quality of oscillations in 
a hypothetical stochastic clock, showing how both mi- 
croscopic oscillations and macrostates are constrained by 



2 



the number of states in the network. The paper con- 
cludes by proposing some experiments which may cast 
direct hght onto the physical realization of these bounds 
on oscillations. 



II. THE LIMITS ON OSCILLATIONS 

To understand the oscillations in the system repre- 
sented by T, we consider the relative contributions of 
different eigenmodes. From the Perron-Frobenius theo- 
rem, all eigenvalues of T have nonpositive real parts, so 
all but the A = equilibrium mode decay away to zero. 
Eigenvectors with nonzero imaginary eigenvalues oscil- 
late in magnitude as they decay. As we see in Eq. |2j 
after a time (Re A^)"^, mode z's contribution to p{t) will 
have substantially diminished. If there is an imaginary 
part to Ai, mode i will oscillate |ImAi/ReAi| times be- 
fore decaying. Because each oscillatory mode will decay 
independent of the other modes, the overall quality of 
oscillations is given by the most oscillatory mode: 
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In a closed system, Q is the upper bound on the to- 
tal number of oscillations arising from an optimal initial 
state of the system. On the other hand, in an open and 
homogeneously driven system, Q describes the coherence 
of those oscillations, in analogy to the quality-factor of 
harmonic oscillators, because external driving may cause 
the system to oscillate indefinitely. This work establishes 
upper-bounds on Q by showing the eigenvalues of T only 
exist in specific regions of the complex plain. 

Karpelevich's Theorem, as clarified by Ito [21] [22], 
states that all possible eigenvalues of an A/'-dimensional 
stochastic matrix with unit spectral radius (|Amax| = 1) 
are contained in a bounded region, R^^ on the complex 
plane, shown in Fig. [l] Rn intersects the unit circle at 
points exp(27rm/6), where a and b are relatively prime 
and < a < b < N. The curve connecting points 
z = e^^^^i/^i and z = e^^^^^/^s described by the para- 
metric equation 

^2 (/I _ LN/biJ = /I LN/6xJ (1 _ lN/b,i ^ 

where s runs over the interval [0,1] and lx/y\ is the 
integer fioor of x/y. For example, the curve that con- 
nects z = 1, corresponding to (ai =0, bi = N)^ with 
z = e^-V^ (a2 = 1, 62 = N) is z{s) = (e^-^^ - l)s + 1. 

The rate matrix from Eq. [l] T, is not a stochastic ma- 
trix. To preserve probability, the sum of each columns 
of T is zero, and the diagonal elements are < 0[23 . To 
transform T into a stochastic matrix, denoted T', divide 
T by the sum of its largest diagonal element and largest 
eigenvalue then add the identity matrix. This transfor- 
mation allows us to write the eigenvalues of in terms 
of the eigenvalues of T: 



a;. = 



max.; IT,; 



max.- I A.- 



1. 
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FIG. 1. The region Rn contains all possible eigenvalues of A/'- 
dimensional stochastic matrices. Region Rn-\-i contains Rn- 
This region is symmetric to the real axis and circumscribed 
by the unit circle. The curves defining each region are given 
by Eq. [4] due to Karpelevich's Theorem. 



Because max^ A^ = 0, and all other eigenvalues of T have 
negative real parts, the most positive eigenvalue of is 
1. Our normalization procedure ensures all other eigen- 
values of are less than 1 and fit within the region Rn 
on the complex plane. Therefore, all of the eigenvalues of 
T fit within the region (max, \Tii \ +maxi \Xi\) - {Rn — 1) . 
Within this region, the maximum number of oscillations 
corresponds to eigenvalues along the line A oc (e±27r2/Ar _ 
1), giving 



Qmax — 



sin(27r/7V) 



cos(27r/Ar) - 1 



= icot(^/7V)<^. (6) 



We can further refine the limit in Eq. [6] using a result 
from Kellogg and Stephens [24 , giving 



Qmax — ^ cot 



< 



eye 



27r ' 



(7) 



where Icye is the longest cycle in the system. 

Up to this point in our proof, we have restricted our- 
selves to systems without any degeneracy in the eigenval- 
ues of T. With degeneracy, as shown in Eq. [2], the time 
dependence of eigenvector j may pick up an extra poly- 
nomial factor, aj(t), with degree less than the degeneracy 
of Aj, which is always less than N — 1. Fortuitous bal- 
ancing of coefficients could allow a p*^-order polynomial 
to add an additional p/2 oscillations. Examining Eq. 2j 
we see that the total maximum oscillation quality can be 



Qmax < ^ + ^ < iV, 



(8) 



where the second term is strictly due to degeneracy ^l] 
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FIG. 2. (a) When the energy landscape has barriers much 
larger than ksT, the system will spend most of its time in the 
minima of the environment. Approximating the continuous 
landscape by discrete states gives the familiar master equation 
kinetics. Here we document two examples of systems with 
unidirectional transition rates. This cyclical system produces 
the maximum Q for any given N. (b) As shown in Eq. [6j 
a system with only two states can not coherently oscillate. 
It produces only random jumps. As the number of states 
in the unidirectional cycle increases (in the same family as 
shown in (a)), oscillations become more coherent and more 
persistent. The spectral density of the unidirectional cycle 
shows a distinct peak which becomes sharper as N increases. 
The transition rates have been normalized by the number of 
states. The Q of the systems are, from bottom to top: 0.75, 
1.46,3.09,6.40,12.7, obtained by fitting Lorentzian functions 
to the peaks. 



III. OSCILLATIONS IN MACRO STATES: 
CHEMICAL CLOCKS 

Oscillations on the discrete state-space are only possi- 
ble when the system in question violates detailed bal- 
ance, which would be the case in a system which is 
driven. Under detailed balance, all eigenvalues of T are 
real [20]. The microstates of this system, which is the 
instantaneous realization of would not show any 

oscillatory dynamics or peaks in the spectral density, 
shown in Fig. [2|a). On the other hand, oscillations in 
macrostates, which are defined as the linear superposi- 
tions {A{t)) = ^iAiPi{t) [26^ (where Ai are some ex- 
perimental observables corresponding to state z), do not 
require the underlying microstates to oscillate. The mi- 
crostate probabilities only need to evolve such that {A{t)) 
oscillates. For example, consider a hypothetical chemical 
clock [2^ described by a cycle of states ^ i28j 



si 



Si 



Si. 



(9) 



If a macroscopic clock advances one tick each time the 
cycle is traversed, adding states to the cycle improves the 
quality of the clock, as shown in Fig[2|b), consistent with 
the limit in Eq. [Sj 

Now consider a more abstract implementation of a 
clock cycle coupled to a fuel reservoir. We define that 
the clock consumes 1 unit of fuel during the transition 
from S£ back to si, so the cycle is an open system driven 



by the fuel reservoir. This accounting method, in ef- 
fect, unrolls the circular cycle in Eq. [9]into a linear chain 
of microstates enumerated by the dyad {/, s^}, where 
/ is the amount of fuel remaining and Si is the state 
of the clock. The macrostate of the clock is {A{t)) = 
Y^ij Aip{f,Si,t) = Y.iAi{p{si,t))f. When / is effec- 
tively infinite, the dynamics of the closed cycle and the 
linear chain are equivalent. As the system moves from 
one state to another, we count time by keeping track of 
the evolution of the macrostate A{t). Therefore, the qual- 
ity of oscillations in the infinite linear chain described by 
the dyad {/, s^} is bounded by Qs < ^, regardless of the 
precise amount of fuel, demonstrating that the oscillatory 
limit of the macrostate is also constrained by the size of 
the network representing the chemical clock in Eq. [9j For 
example, take = {1 if mod2i = 0; otherwise}, mean- 
ing i = 2. The evolution of {A{t)) is shown in Fig. [sj^a). 
The total number of states, does not effect the quality 
of the clock. However, if we change Ai to 

A^ = {1 if modpi = 0; otherwise}. 

Fig. [sjb) shows increasing oscillation in {A{t)) with in- 
creasing p. That is, adding more states to the clock cir- 
cuit directly increases its accuracy. 

Texts exploring chemical oscillations state that nonlin- 
earity is a requirement for oscillations. In fact, nonlin- 
earity is a shorthand for describing extremely large sys- 
tems [13 . Under conditions of detailed balance, systems 
must consume some sort of fuel to sustain oscillations. 
If we consider the fuel-free states as being an abstract 
engine with N states, the combined engine-fuel system is 
describe by being in one of N different states and having 
/ units of fuel remaining. Therefore, a fuel reservoir can 
allow a total number of oscillations ~ fN/27T. Because 
we can approximate an extended fuel-engine network as 
the engine network alone combined with rates that ac- 
count for the chemical potential of the fuel, Eq. |8] implies 
that the number of inherently unique states of the engine, 
absent fuel consumption, will constrain the possible reg- 
ularity of reciprocal motion [9 . 

Similarly, if the system is driven by oscillations in mul- 
tiple parameters or species, we can again parameterize 
the state of the system based on the population of each 
component. However, most chemical dynamics are mod- 
eled using continuous variables, not discrete numbers of 
states. The microscopic description of the system, com- 
prised of a discrete number of states, is connected to 
the continuous mass-action approximation of chemical 
dynamics by a system size expansion described by Van 
Kampen [13]. Take, as an example, the multidimen- 
sional oscillating chemical reaction called the Brussela- 
tor. By expanding the mass-action Brusselator into a 
discrete state-space, the system size expansion parame- 
ter determines the length of the largest cycle [29l [30] . 
Indeed, multiple authors have observed that the qual- 
ity of the Brusselator limit cycle scales with system size, 
consistent with this work [28-33 . For example, Gaspard 
et al. derived that the diffusion in phase-space of the 
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frequency (arb. units) frequency (arb. units) 

FIG. 3. Although the hnear set of states given by Eq. [9]does not have any imaginary eigenvalues, macrostates can oscillate. 
Macrostates are defined as {A) = Y.i ^^^P^{-^^)^ In this case, = {1 if modpi = 0; otherwise}, (a) Dynamics for different 
values of N for fixed p = 2 (inset) Identical dynamics, but with scaled rates so expected traversal times are the same, (b) 
Increasing p increases the number and coherence of oscillations given fixed N — 100, demonstrating the limit in Eq. [7| 



Brusselator, which corresponds to amount of damping in 
the macroscopic observables, is inversely proportional to 
the size of the system [3^ [33] . This observation is not 
merely coincidence, but a fundamental efficiency limit of 
the master equation. 



IV. CONCLUSION 

The bounds on oscillations can play a key role in inter- 
preting experimental observations by determining a min- 
imum number of underlying states. For example, the os- 
cillation of fluorescence wavelength in fluorescent protein 
GFPmut2 remains unexplained |34H36] . After applica- 
tion of a denaturant, the ionic state of the fiuorophore can 
switch up to Q ~ 50 times with high regularity, observed 
as oscillations in the emission wavelength [34J. Because 
Eq. [T] bounds the number of states involved in the oscil- 
lation to be at least 27r times larger than Q, this predicts 
that the oscillations are driven by large-scale rearrange- 
ment of the numerous hydrogen bonds in the /3-barrel, 
not merely exchange between the few amino acids di- 
rectly connected to the fiuorophore. If the protein were 
to be mutated to alter the number of bonds in the 
barrel, we predict that we should see a corresponding 
alteration in the number and quality of observed oscilla- 
tions. 

The quality bounds proved hear are universal. Be- 
cause the master equation is used in nearly every branch 
of science, the dynamics being modeled need not be phys- 
ical |2^. For example, it could be money held by a 
bank [^, packets of data on the internet [37 , agents 
traversing a network [38 , or the populations in an ecosys- 
tem [39]. The oscillation limit could also be probed 



experimentally with sculpted landscapes using optical 
tweezers [40 . As a probe bead jumps from trap to trap, 
the energy landscape in unoccupied traps is sculpted to 
simulate an arbitrarily large designer network of discrete 
states. Ultimately, the current results are a fresh ap- 
proach to analyzing the dynamics of discrete systems, 
and it serves as a new design principle for those seeking 
to engineer oscillations. 

This efficiency limit of oscillations has obvious impli- 
cations on how well a high-dimensional system can be 
numerically approximated by a smaller system. The ap- 
proximation will only be successful if the relevant eigen- 
values of the larger system lie within the allowed region 
of the smaller system. However, the inverse stochastic 
eigenvalue problem has not yet been solved, so we can 
not know a priori if a stochastic matrix exists for a given 
set of eigenvalues, even if they all reside within the al- 
lowed region [41 . This fact prevents us from constructing 
the opposite bounds, the conditions for a minimum num- 
ber of oscillations. Hopefully, future results will further 
constrain the present bounds, and we may gain deeper 
insight into the necessary conditions for creating oscilla- 
tions. 
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